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Abstract 

We analyze an effective model on a square lattice with two types of forward scattering interac- 
tions, which, respectively, drive ferromagnetism (FM) and electronic nematic order via a d-wave 
Pomeranchuk instability (dPI). The FM and dPI in general compete with each other and they 
are typically separated by a first order phase boundary in the plane of the chemical potential and 
temperature. Nevertheless there is a parameter region where the dPI occurs inside the FM phase, 
leading to their coexistence. We also study the effect of a magnetic field by choosing a chemical 
potential where the ground state is paramagnetic without a field. In this case, instead of FM, the 
dPI competes with a metamagnetic instability. The latter occurs above a threshold strength of 
the FM interaction and otherwise the dPI is stabilized with a dome-shaped phase diagram in the 
plane of a magnetic field and temperature. The FM interaction shifts the center of the dome to a 
lower field, accompanied by a substantial reduction of the field range where the dPI is stabilized 
and by an extension of the first order part of the transition line, although the maximal critical 
temperature does not change. The experimental phase diagram of the bilayer ruthenate Sr3Ru207 
can be well captured by the present theory. 

PACS numbers: 71.27.-l-a, 71.18.+y, 75.25.Dk, 74.70.Pq 
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I. INTRODUCTION 



In the nematic liquid crystal,- rodlike molecules have a preferred orientation. This state is 
characterized by breaking of orientational symmetry, retaining other symmetries of the sys- 
tem. Electronic analogues of the nematic liquid crystal attract much interest. Electrons have 
spin and the direction is defined in spin space. Using spin degrees of freedom, a spin nematic 
state is studied in quantum spin systems.- 1 ^ Electrons also have orbital degrees of freedom. 
With orbital order such as an occupation difference between the d yz - and <i 2:r -orbital in a 
(i-electron system, electrons may break orientational symmetry without any additional sym- 
metry breaking, leading to an orbital nematic state.- 1 ^ Ferropnictides are possible materials 
for such a state.— On the other hand, the orientation cannot be defined for charge itself. 
However, a nematic state can be realized by using a charge degree of freedom. Two routes to- 
ward a charge nematic state are proposed. When the system is close to a charge stripe order, 
namely one- dimensional- charge order, where both translational and orientational symmetry 
are broken, fluctuations of charge stripes may restore the former but the latter may be still 
broken.- The charge nematic order can be obtained also without invoking charge stripes. 
It was found theoretically that the two-dimensional t-J— and Hubbard 10 models exhibit a 
tendency toward a d-wave Pomeranchuk— instability (dPI). In this state, the Fermi surface 
expands along the k x direction and shrinks along the k y direction, or vice versa, whereas in 
a real space representation the nearest neighbor hopping is effectively enhanced along one 
direction and suppressed along the other direction. 

The <£PI was extensively studied not only in the t-J—^~— and HubbarcU^£r— models, 
but also in phenomenological models,— >2i a model with central forces,— >22 general Fermi 
liquid schemes,— 1 ^ and continuum (not lattice) models.— ~— Mean-field theory of the (iPI— i2i 
showed that the dPI occurs around van Hove filling with a dome-shaped transition line. 
Typically the transition is second order at high temperature and changes to first order at 
lower temperature. The end points of the second order line are tricritical points. The mean- 
field phase diagram is characterized by a single energy scale, similar to the BCS theory of 
superconductivity, and thus various universal numbers were found.™ 

Fluctuations of the <iPI suppress the first order transition obtained in mean-field theory 
and when they are strong enough, the transition changes to be continuous even at zero 
temperature, leading to a quantum critical point.— ^ At the quantum critical point, dPI 
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fluctuations lead to a non-Fermi liquid ground state.— ^ At finite temperatures close to the 
dPl, thermal fluctuations become dominant. They turned out to truncate the original Fermi 
surface, leading to a Fermi-arc-like feature.— 

Signatures of nematicity were observed in cuprate superconductors. Neutron scattering 
measurements revealed a strong anisotropy of magnetic excitations in momentum space.— - — 
The anisotropy showed strong temperature and doping dependences, which are well cap- 
tured in terms of the competition of the dPl and singlet pairing formation.— & Transport 
measurements also revealed a very strong anisotropy of the Nernst coefficient, 4 - which was 
interpreted ClS db SI gnature of charge nematic order.— 

There is growing evidence that the bilayer ruthenate Sr 3 Ru207 (Sr327) exhibits a dPl 
in a strong magnetic field.— - — In fact, many features observed in experiments were well 
understood in terms of the dPl, for example, the metamagnetic transition,— the enhancement 
of the residual resistivity— the bilayer effect,— 1 ^ the suppression of the critical temperature 
by impurities,— and the spin-orbit effect.— Furthermore, the experimental phase diagram 
is very similar to that obtained in mean-field theory— In particular, it was found that the 
mean-field phase diagram is characterized by a single energy scale even in the presence of a 
magnetic field.— Therefore there exist various universal ratios for a given chemical potential, 
which can be compared directly with experimental data. Although several universal ratios 
agree with the experimental data, ratios of the characteristic temperature and field give one 
order of magnitude smaller than the experimental ones.— 

This apparent inconsistency cannot be resolved by invoking different choices of param- 
eters. The key may lie in the set of experimental indications that Sr327 is located close 
to a ferromagnetic instability: a large Wilson ratio,— a uniaxial-pressure-induced ferromag- 
netic transition, 56 - and the presence of ferromagnetic fluctuations observed by the inelas- 
tic neutron scattering,— the nuclear spin-lattice relaxation rate,— and thermal expansion 
measurements.— Moreover several band calculations^ 1 ^ for Sr327 (without a field) suggested 
that the system is close to ferromagnetism (FM). Hence the presence of a ferromagnetic in- 
teraction is quite plausible in Sr327. In fact, early theoretical work— ^ for Sr327 focused on 
the role of ferromagnetic interactions, especially in the context of a metamagnetic transition 
observed in experiments.— 

In this paper, we develop a mean-field theory by taking two types of forward scattering 
interactions, which drive the <fPI and FM, respectively, into account. In the context of Sr327, 
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it is interesting to explore how the mean-field phase diagram of the dPI obtained previously 
is modified by the presence of a ferromagnetic interaction and how well the experimental 
phase diagram of Sr327 is captured. Furthermore, the interplay of the dPI and FM is 
interesting in its own right. While FM is an instability in the spin channel whereas the <iPI 
is in the charge channel, both are instabilities in the particle-hole channel of q=0 and do 
not break translational symmetry. Several theoretical analyses of microscopic models^ 1 ^ 1 ^ 
actually suggested the presence of a ferromagnetic instability, which competes with the dPI. 
Therefore in a more general setting we study the interplay of the <iPI and FM, and clarify 
possible scenarios in such a coupled system. 

We propose an effective model, suitable to address the interplay of the dPI and FM, and 
derive resulting phase diagrams. In Sec. II, we introduce a forward scattering model and 
present results in Sec. Ill by separating two cases: i) zero magnetic field (h = 0) and ii) 
finite magnetic field (h ^ 0). The latter case is relevant to Sr327. Conclusions follow in 
Sec. IV. 

II. MODEL 

To analyze a coupled system of the <iPI and FM, we consider the following Hamiltonian 
on a square lattice, 

H = Ho + Ti^ + Tim + Hz ■ (1) 
The first term "Ho is the kinetic term, 

^o = J>£-/i)i CT c kCT , (2) 

kcr 

where c kCT (cko-) is a creation (annihilation) operator of an electron with spin a and momen- 
tum k; fi is the chemical potential. The electron dispersion is given by 

e k = — 2t(cos k x + cos k y ) — At' cos k x cos k y (3) 

with t and t' being the nearest and second nearest neighbor hopping amplitudes, respectively. 
The second term "H^ is a forward scattering interaction driving a dPI, 

n ^ = ~2N ^ d kd k '4 CT c kCT 4, CT ,c k v , (4) 

kkW 
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where the coupling constant is positive, is a <i-wave form factor such as = cos k x — 
cos k y , and N is the total number of lattice sites. This term describes the ci-wave weighted 
density-density interaction with zero momentum transfer, which was obtained in microscopic 
models such as the t-J— and Hubbard^ 1 ^ models. 

The third term % m describes an Ising ferromagnetic interaction, 




where g m {> 0) is a coupling constant and a = +1 and —1 for up-spin and down-spin, 
respectively. This interaction is obtained by focusing on the spin-spin interaction with a 
spin quantization axis parallel to the z direction and by extracting a scattering process with 
zero momentum transfer. Therefore the interaction described by 1-L m is appropriate when 
the system has a strong spin anisotropy as well as dominant forward scattering processes 
of electrons. The interaction of l-i m is also obtained by considering a mean-field analysis 
of spin rotational invariant interactions. For instance, in the case of the Hubbard onsite 
interaction U J2i n «t n 4> our coupling constant is given by g m = 2U. 
The last term Hz is the Zeeman energy, 

Hz = -2^2 (T( Lt c ^- ( 6 ) 

ko 

Here h is an effective magnetic field given by h = jj/ieif , with q being a g factor, fi B is the 
Bohr magneton, and if is a magnetic field. 

The terms of and T-L m describe pure forward scattering interactions of electrons. Thus 
fluctuations around the mean-field vanish in the thermodynamic limit. In other words, 
mean-field theory solves our Hamiltonian exactly in the limit of iV — > oo. 

The order parameter of the dPl is defined by 

kcr 

This quantity becomes finite only if the system breaks square lattice symmetry because of 
the presence of the (i-wave form factor. FM order is defined by 

m = f^5Z a ( c ko c ko), (8) 

kcr 

where we include the coupling constant g m in the definition of m; while the magnetization 
is then given by m/g m , we may refer to m as magnetization, as long as no confusion occurs. 
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We decouple the interaction terms (@| and (jSJ) by introducing the order parameters and 
m, and obtain the mean-field Hamiltonian, 

U M f = ^ai CT c k(J + + -^-(f 1 , (9) 

kcr ^ m ^ 

where the renormalized dispersion is given by 

Cka = e° - |(m + /i) - 40 - • (10) 
The grand canonical potential per site at temperature T is obtained as 

w = ~ X) Ml + e' WT ) + + ^ • (11) 

N ^ 2g m 2g (j) 

The stationary condition of u with respect to (f) and m leads to the self-consistent equations 

^ = ^E^/(w, (i 2 ) 

kff 
kcr 

which we solve numerically. Here /(£k<r) is the Fermi function. 
III. RESULTS 

We fix gr^/t = 1 throughout this paper unless otherwise noted and explore how the phase 
diagram of the riPI changes with increasing the FM interaction g m . We first study the case 
of h = and then that of h ^ 0. As a band parameter, we choose t'/t = 0.35, which was 
used for the study of Sr327.— ^ Since the presence of t' turns out to play a crucial role to 
understand phase diagrams for h = 0, we also study the case of t' — for h = 0. Hereafter 
we set t = 1 and all quantities with dimension of energy are in units of t. 

A. Results for h = 

1. Evolution of phase diagrams with increasing FM interaction 

Figure [T] shows a sequence of phase diagrams for g m < 7.0 in the plane of the chemical 
potential /x and temperature T. Because of the competition with the dPl, no FM instability 
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occurs at least up to g m = 6.0 [Fig. [T] (a)] and the phase diagram is occupied only by the 
dPl. As already clarified previously,—* 2 ^ the dPl occurs below a dome-shaped transition line, 
with a maximal T c near the van Hove energy (/i v h = At' = 1.4); a deviation from /i v h is 
due to the presence of t', which breaks particle-hole symmetry. The phase transition is of 
second order at high temperature (T 2nd ) and of first order at low temperature (T c lst ). The 
end points of the second order line are tricritical points (T* 
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FIG. 1: Phase diagram in the (/U,T) plane for a sequence of couplings g m . Transition from the 
paramagnetic to ordered phase is a second order (T c 2nd ) at high T and a first order (T c lst ) at low 
T; T* n is the temperature at a tricritical point. A dashed line (Tl^) denotes the first order phase 
boundary between the <£PI and FM, which appears in (b) and (c). 



For g m = 6.5, the FM interaction becomes strong enough to realize FM near the edge on 
the side of a high chemical potential [Fig. [1] (b)]. The transition from the paramagnetic to 
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FM phase is second order at high temperature, but the second order line ends at a tricritical 
point and changes to a first order line at low temperature. This feature is the same as the 
transition between the paramagnetic and dPl phase. The boundary of the dPl and FM is 
characterized by a first order transition (T^). 

As shown in Fig. [T] (c), this first order phase boundary shifts to the middle of the phase 
diagram for g m = 7.0 and the FM becomes more stable. The order parameters are plotted 
as a function of fi in Figs. [2] (a) and (b) at T = 0.01 and 0.20, respectively At a low 
temperature (T = 0.01), and m show a jump at \i ~ 1.05 and 1.81, respectively, because 
of a first order transition from the paramagnetic phase. The <iPI changes to the FM via a 
first order transition at fi 1.45 and there is no mixing of and m. At a high temperature 
(T = 0.20), on the other hand, (ft and m develop continuously at fi « 1.10 and 1.72, 
respectively. The transition between the dPl and FM is however still of first order. 
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FIG. 2: (Color online) \i dependence of 4> and m at T = 0.01 (a) and 0.20 (b) for g m = 7.0. 



As expected, with further increasing g m , the first order boundary between the <iPI and 
FM shifts to a lower chemical potential. In fact, as shown in Fig. [3] (a), the c£PI is realized 
only near the edge of the dome for g m = 7.8. However, qualitative changes occur in the 
phase diagram. First, the coexistence of the dPl and FM is stabilized inside the FM phase 
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near the edge of the first order line of the FM around fi = 2.04. This region is magnified in 
Fig. [3] (b). The transition from the FM to the coexistence is first order at low temperature 
and becomes second order at high temperature. While one end point of the second order line 
at fi w 2.037 is a tricritical point, the other end point at \i ~ 2.045 is just a point touching 
the first order line of the FM. There is a direct first order transition from the paramagnetic 
phase to the coexistence around = 2.05. Second, an additional FM phase appears in 
2.52 < < 2.6 as shown in Figs. [3] (a) and (c). This FM comes from the enhancement of 
the density of state at the band edge of /i = 2.6. A first order transition occurs only on the 
side of a lower chemical potential and the second order line disappears at the band edge. 
This band-edge FM is realized for 7.6 < g m < 8.0. 




FIG. 3: Phase diagram in the (fi,T) plane for g m = 7.8. The regions near \i = 2.04 and 2.55 are 
magnified in (b) and (c), respectively. 

For g m = 8, as shown in Fig. H] (a), the FM becomes dominant and a pure dPl phase is 
not stabilized. Instead the <iPI is realized in coexistence with the FM around fi = 2.06, as 
magnified in Fig. @] (b). In contrast to the case of g m = 7.8 [Fig. [3] (b)], the phase boundary 
of the coexistence is well separated from the first order line of the FM, leading to a phase 
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(a) 




FIG. 4: (a) Phase diagram in the ((J*,T) plane for g m = 8.0. The region of the coexistence around 
fjL = 2.06 is magnified in (b). (c) [i dependence of m at T = 0.001. Two successive jumps around 
fi = 2.06 are magnified in the inset. 

diagram very similar to that of the pure dPl [Fig. [1] (a)], but with a significant extension 
of the first order portion of the transition line; the reason for this will be explained later in 
terms of Eq. ( 12Tj) . The magnetization m is plotted as a function of fi in Fig. H] (c) at low 
temperature. After the first order FM transition at fi ~ 1.05, the value of m increases with 
increasing /i and forms a cusp at fi ~ 1.45 where the density of states of up-spin electrons 
is fully occupied and the system changes to a half-metallic state. For fi > 1.45, m decreases 
since electrons with down-spin increase whereas the electron density of up-spin remains 
unity. At fi w 2.05 and 2.07, m exhibits a jump [see the inset of Fig. H] (c)] because of the 
presence of the coexistence of the <iPI and FM, which occurs via a first order transition at 
low T. The magnetization m vanishes discontinuously at fi ~ 2.15, but appears again with 
a jump at /i ~ 2.42 because of a first order transition associated with the band-edge FM. 



10 



The magnitude of m decreases monotonically and vanishes at the band edge of /i = 2.6. The 
system becomes a band insulator for \i > 2.6. 

With further increasing g m (Fig. [SJ, the band-edge FM is absorbed into the main FM 
phase. A first order phase transition then occurs only on the lower side of /i. Inside the 
FM, the coexistence of the <iPI and FM is stabilized up to g m = 9.8. Figure (a) is 
the representative phase diagram computed for g m = 9. In Fig. [5] (b) the region of the 
coexistence of the dPl and FM is magnified. This phase diagram is very similar to that for 
g m = 8 [Fig. H] (b)] with the same maximal T c , but with a further extension of the first order 
transition line. For g m > 9.8, however, the coexistence is replaced by a first order transition 
associated with a jump of the magnetization, namely a metamagnetic transition inside the 
FM, as denoted by solid squares in Fig. \5\ (c). The magnetization is potted as a function of 
H at low T in Fig. [5] (d). The jump at fi « 2.23 comes from the metamagnetic transition. 
The cusp at /i ~ 0.81 indicates that the up-spin band is fully occupied in /i > 0.81, where 
the system becomes half-metallic. 
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FIG. 5: (a) Phase diagram in the (//, T) plane for g m = 9.0. The region of the coexistence is 
magnified in (b). (c) Phase diagram for g m = 10. T m denotes the position where a metamagnetic 
transition occurs, (d) \i dependence of m at T = 0.001 for g m = 10. 
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2. Discussions 



The coexistence of the dPl and FM is stabilized even for g m ^> (Figs. [3H5]). This is 
because of the presence of the van Hove singularity. After performing explicit calculations 
up to g m = 10, we confirm the van Hove singularity due to the down-spin band (m > is 
assumed) inside the FM phase for g m > 7.8. Around the van Hove filling, either the dPI or 
a metamagnetic transition occurs in our model, depending on energetics. We find that the 
coexistence of the dPI and FM is more favorable for 7.8 < g m < 9.8 and the metamagnetic 
transition for g m > 9.8. 




FIG. 6: Phase diagram in the (/x,T) plane for a sequence of couplings g m by setting t' = 0. The 
phase diagram is occupied by the dPI in g m < 8.84 (a) and by the FM in g m > 8.87 (d). In a tiny 
range of g m [(b) and (c)], both FM and dPI are realized, but separated from each other by a first 
order boundary; the line of Tj^ appears only in (b) and (c). 

Our results shown in Figs. HHS] are very asymmetric with respect to the van Hove energy 
of the bare dispersion, which is given by /i v n = 4i' = 1.4. This is because the presence of 
t' breaks particle-hole symmetry. In fact, for t' = 0, the phase diagram becomes symmetric 
with respect to the axis of /x = 0. For < g m < 8.84, the dPI is stabilized and no FM 
is realized [Fig. [6] (a)]. For g m > 8.85, however, the dPI starts to be replaced by the FM 
phase from a higher temperature [Fig. |6] (b)] and is stabilized only around /i — at low T 
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for g m = 8.86 [Fig. |6] (c)]. The dPl disappears already for g m = 8.87. The change from the 
dPl [Fig. |6] (a)] to the FM phase [Fig. [6] (d)] occurs in a very small range of g m . In contrast 
to the case of Figs. El EH and El no coexistence of the dPl and FM is stabilized. Furthermore 
a band-edge FM does not appear. 

Our results for h = are summarized as follows: i) in < g m < g m x, only the dPl phase 
is realized, ii) in g m i < g m < g m 3, both c£PI and FM are stabilized, but they are separated 
from each other by a first order transition line, iii) in g m 2 < g m < 9mA, t ne coexistence with 
dPl occurs inside the FM phase, and iv) in g m ^ < g m , only the FM is realized. We have 
obtained g m \ ~ 6.5, g m2 ~ 7.8, g m3 ~ 7.9, and g mi ~ 9.8 for t' = 0.35, leading to rich phase 
diagrams as shown in Figs. [TJ El IU and El For t' = 0, on the other hand, we have obtained 
9mi ~ 8.84, g m2 = g m 3 = g m ± ~ 8.87. As a result, a phase diagram is occupied by either the 
dP\ or FM except for a tiny range of g m . 



B. Results for h ^ 

Next we examine the effect of a magnetic field, motivated by the experimental indication 
that Sr327 is paramagnetic in zero field and exhibits a nematic instability around 8 Tesla.— _ — 
Fixing the chemical potential /i = 1 and taking the field as a tuning parameter, we study 
how the phase diagram of the <iPI evolves with increasing the ferromagnetic interaction. 

Figure [7] (a) is a set of phase diagrams of the dPl in the plane of a magnetic field and 
temperature for a sequence of g m , showing four characteristic features: with increasing g m , 
i) the dPl occurs in a lower field, ii) the field range where the dPl is stabilized shrinks 
substantially, iii) the first order part of the transition line extends and tricritical points are 
pushed up to higher temperatures, but iv) the maximal T c does not change. 

To understand these features, we consider a magnetic field h v u, at which the cr-spin band 
touches the van Hove energy, and the <iPI is expected around that. From Eq. ( flOl) . /i v h 
fulfills for = the relation, 

a(m + h vH ) 

+ /i = /ivH , (14) 

and the corresponding relation for the other spin band should be —aim + /i v h)/2 + /x = 
2/i — /i v H, where /i v n = 4t'. Since /i is fixed in our case, we obtain 



2-vH 



2|/i — /i vH | — m . (15) 
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FIG. 7: (Color online) (a) Phase diagram in the (h,T) plane for a sequence of couplings g m ; the 
value of g m is indicated near the maximal T c . The dPI is stabilized inside the dome for each g m . 
(b) h dependence of (ft at T = 0.001 for a sequence of g m . (c) h dependence of m for g m = 7 at 
T = 0.001. The corresponding result for g^ = is also plotted. 



While the magnetization is not fully linear in field in the entire field range we consider, we 
may invoke the equation obtained in linear response theory, 

m/g m ~xhvn, (16) 
= 7^^, (17) 

1 - 9mX0 

where \ is the full magnetic susceptibility, which is expressed by the non-interacting magnetic 
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susceptibility xo as shown in the second line; the presence of g m on the left-hand side is due 
to our definition of m [see Eq. (jSJ)]. We then obtain 

h vll = 2(1 - g m Xo)\n - AMil , (18) 

that it, the value of h v n is reduced with increasing g m . Since the dPl occurs around the van 
Hove energy, the dPl should be realized around a lower field with increasing g m . 

Equation (fl6|) is a rough approximation especially near a metamagnetic transition and 
the resulting Eq. (ITS]) should be taken as qualitative understanding. To get quantitative 
understanding, we solve Eq. ( TT3]) numerically under the condition of = and Eq. (IT4]) . 
We then obtain /i vH ~ 0.61,0.41,0.22,0.12 for g m = 2,4,6,7, respectively; for g m = 0, on 
the other hand, h v n = 2|/z vH — fJ>\ — 0.8 since m = 0. The dPl indeed occurs around those 
fields in Fig. [3 

The range of a magnetic field where the dPl is stabilized becomes narrower for a larger 
g m . As seen in Eq. ( ITU]) , the sum of m and h plays a role as an effective field. Since m 
becomes more susceptible to a field as g m becomes larger and furthermore m is proportional 
to g m in our definition [Eq. (JSJ], the value of h to stabilize the dPl is necessarily reduced. 

The first order transition line extends with increasing g m . To understand this, we expand 
the free energy Eq. (1111) with respect to the order parameter of the dPl around = 0, 

oj{4>\ m) — u)(0; m) = ^a 2 2 + -jy a 40 4 H • (19) 

The coefficients of a 2 and a 4 are obtained as 

^ V ka / 

1 r4 J-W/fO \ o„ (2IV S^ko- °"^k/ (^kg)) / Q1 x 

°4 - T7 2^ d ki l4k<xJ - 6 9m - 5m v 77770 \ ' 

where ^ ko . = e k — °"( m 2 +fc ) — ^ and /', /", /'" are the first, second, third derivative of the Fermi 
function. When a 4 becomes negative, a first order transition can occur. The second term 
on the right-hand side of Eq. (12ip originates from the dependence of m. The denominator 
of this term is positive close to the <iPI and the numerator becomes in general finite when 
the spin symmetry is broken. Hence the second term is negative for h ^ 0. Furthermore the 
second term is proportional to g m . Therefore the presence of the second term in Eq. (I2TT) leads 
to an extension of the first order transition line of the dPl and this effect becomes stronger 
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for a larger g m . The same argument explains the extension of the first order portion of the 
transition line in Figs. H] (b) and [5] (b), since the second term of Eq. ( 12 lj) becomes negative 
also in the FM phase. 

A second order transition is given by the condition a 2 = 0. Since /1 is fixed, the quadratic 
term a 2 is a function of h = m + h. Suppose the maximal T c is obtained at h m3X , there 
can exist a field h and a magnetization m, which give the same value of h max for a different 
g m , although the values of m and h themselves depend on g m . This is actually the case up 
to g m = 7.8, leading to the same maximal T c in Fig. [7] (a). A similar consideration also 
explains the same maximal T c in Figs. S] (b) and (b). Keeping in mind that our system is 
half-metallic in the range of \i where the coexistence is stabilized [see the discussion about 
Fig. IH (c)] and thus only the down-spin band is active, the coefficient a 2 becomes a function 
of the quantity ft = = ^ + fi for h = 0. We confirm the same value of ft at the maximal T c 
in Figs. H] (b) and [5] (b), respectively, which necessarily yields the same maximal T c . 

In Fig. [7] (b), the order parameter of the dPl is plotted as a function of h for a sequence 
of g m at low T. Because of two first order transitions at low T [Fig. [7] (a)], the order 
parameter exhibits two jumps. Interestingly the maximal value of <fi does not depend on g m . 
This feature is easily understood from Eqs. f lTUj) and f fT2|) . The right-hand side of Eq. f fT2|) 
depends on the quantity ftk a = |(m + h) + d^cp for a fixed /i. Suppose the maximal value of 
0, say max , is obtained at h = /i max for g m = 0, namely for m = 0. Even when g m is turned 
on, the same value of max is obtained as long as m and h fulfills the equation 

m + h = /i max . (22) 

This equation may hold unless the value of m becomes as large as /i max . We can check that 
Eq. (|2"2"|) indeed holds up to g m ~ 7, leading to the same maximal value of (ft for g m = — 7. 

In Fig. [7] (c), the magnetization is plotted as a function of h at low T. Because of first 
order transitions at low T, the magnetization exhibits two successive jumps. It is instructive 
to recognize that there could occur a metamagnetic transition at h w 0.12 if the coupling g^ 
would be turned off, indicating the underlying competition of the dPl and a metamagnetic 
transition. We can check that the dPl overcomes the metamagnetic transition up to g m = 7.9. 

For g m > 8, on the other hand, the metamagnetic transition becomes dominant and 
the magnetization exhibits a single jump as shown in Fig. [8] (a). The Landau free energy 
is plotted in Fig. (b) as a function of m at h — 0.023, just below the metamagnetic 
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FIG. 8: (Color online) (a) h dependence of m for g m = 8 at T = 0.001. (b) Free energy as a 
function of m at h = 0.023 and T = 0.001 for g m = 8. The value of eft which minimizes the free 
energy at each m is also plotted. 



transition; the order parameter is optimized to minimize the free energy at each m. 
There are three local minima. Two local minima, where = is stabilized, are associated 
with the metamagnetic phenomenon. The other local minimum, at which becomes finite, 
corresponds to a solution of the dPl. This solution, however, does not give the absolute 
minimum and thus the dPl does not occur. 

When g m exceeds 8.25, the FM occurs even for h = 0. In this case, neither a metamagnetic 
transition nor a dPl occurs by applying a magnetic field. 

The effect of a ferromagnetic interaction on the <iPI for h ^ can be summarized as 
follows: i) the dPl occurs in a lower magnetic field, ii) the field range where the <iPI is 
stabilized becomes narrower, iii) the first order part of the transition line extends, and iv) 
the ciPI and a metamagnetic transition compete with each other and the former is realized 
up to g m ~ 8, and the latter for 8 < g m < 8.25 for the present choice of parameters. 
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IV. CONCLUSIONS 



We have studied a two-dimensional electron system, where electrons interact with each 
other via interactions favoring a dPl and FM. In the absence of a magnetic field, we have 
obtained rich phase diagrams. The dPl and FM typically compete with each other. In 
fact, while both <iPI and FM can be realized simultaneously, they are separated by a first 
order phase boundary. Nevertheless it is possible that the dPl is stabilized inside the FM 
phase, leading to their coexistence. The presence of t', leading to a breaking of particle-hole 
symmetry, plays an important role. For t' = 0, either the <iPI or the FM is typically realized 
in the plane of the chemical potential and temperature, and coexistence is not stabilized. 
We have also studied the effect of a magnetic field, motivated by the experimental indication 
that Sr327 is in the normal state without a magnetic field and exhibits a nematic instability 
by applying a field. In this case, instead of FM, the <fPI competes with a metamagnetic 
transition. The latter occurs above a threshold strength of the FM interaction and otherwise 
the dPl is stabilized with a dome-shaped transition line around the van Hove energy in the 
plane of a field and temperature. With increasing the FM interaction, the center of the dome 
shifts to a lower field, accompanied by a substantial reduction of the field range where the 
dPl is stabilized and by an extension of the first order part of the transition line, although 
the maximal T c does not change. 

It might seem that the interaction strength of g m is considered up to a too large value 
(g m ~ 10) in our study. However, this seemingly large value is due to our definition of g m 
in Eq. (J5J) where a factor of (1/2) 2 originating from spin is not absorbed into the definition 
of g m . 

A typical feature of the dPl is that its mean-field phase diagram is characterized by 
universal ratios.— 1 ^ In the model solved in Ref. 



54 



several universal ratios reasonably agree 



with experimental values, but ratios of temperature and a magnetic field come out one order 
of magnitude smaller than the experimental data. For example, in experiments, T^ n /h tT i ~ 
0.6fc_B/(0.15g;U R ) = 6g _1 « 3 if g = 2, whereas theoretically we obtain T c tri //i tri ~ 0.3 for 
9m — (Ref. 1671 ); here h tri is the field at a tricritical point measured from the van Hove 
energy. However in the presence of a ferromagnetic interaction, we have found that only the 
scale of a magnetic field is substantially reduced while the temperature scale is not. As a 
result, from Fig. [7] (a), we obtain T c tn //i tri ~ 2 and 7 (Ref. |67J) for g m = 6 and 7, respectively. 
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The ratio of T* ri //i tr i is substantially modified by a FM interaction to become comparable 
to the experimental one. The large value of g m indicates that the system is close to the FM 
instability for h = 0, the same situation as in Sr327.— ~— 

The FM interaction pushes up T c tri to a higher temperature, but the maximal T c does 
not change. As a result, other ratios such as T* n /T c vH , where T c vH is T c at the van Hove 
energy, now becomes slightly larger than the experimental value, although it showed better 
agreement with experimental data in the model with g m = 0. 54 However, this may be easily 
improved by invoking weak fluctuations associated with the dPl, since it was shown^ 2 -^ 
that fluctuations suppress T* n stronger than T^ n . Therefore the ratios in the experimental 
phase diagram of the c£PI are well understood by the presence of a FM interaction tuning 
the system close to the FM instability, and by weak dPl fluctuations. 

The lines of first order phase transitions tilt outward in the experimental phase diagram,— 
indicating that the entropy inside the <iPI phase is larger than that in the normal state.— 
This counterintuitive phenomenon is not captured in the present theory. This inconsistency 
may be explored further in terms of the interplay of ferromagnetic fluctuations and the <iPI 
by going beyond the mean-field model. 

While we have analyzed a single band model, Sr327 is a t 2g system and orbital nematic 
order may provide another possible scenario.-^ Since there are interactions among different 
orbitals, the <iPI is expected to generate orbital nematic order, or vice versa. It is an open 
question which is the driving force for nematicity observed in Sr327. 
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